Effect of van Hove Singularities on a Spin Liquid 
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We determine the properties and leading instabilities of a spin liquid with a Fermi surface passing 
near a van Hove singularity. Our study is motivated by recent photoemission experiments on high 
Tc cuprates in which it is found that for the optimally doped material the experimental Fermi 
surface passes near a van Hove singularity, while for underdoped materials, a pseudogap in the 
electron spectral function is formed in the vicinity of the van Hove point. We show theoretically 
that proximity to the van Hove singularity suppresses the inelastic scattering due to the gauge 
field and permits the formation of a d-wave RVB state in which the gap exists only near the van 
Hove points while finite regions of the Fermi surface remain gapless. This d-wave pairing provides 
a natural explanation of the pseudogap observed in photoemission. We also discuss the relation 
of the pseudogap observed in the spectral function to the pseudogaps observed in the magnetic 
susceptibility. 



I. INTRODUCTION 

In this note we report results of a theoretical study of a 
"spin liquid" with a van Hove singularity near the Fermi 
surface. By the term "spin liquid" we mean a liquid of 
charge spin 1/2 fermions filling a large Ferrni sea and 
coupled by a singular gauge field interactionldu and by 
additional non-singular interactions. We determine ex- 
actly the effect of the van Hove singularity on the fermion 
gauge-field physics and treat the additional interactions 
by a "leading logarithm" renormalization group analysis. 

The problem of a spin liquid with a van Hove singular- 
ity is of interest on experimental and theoretical grounds. 
The experimental motivation for describing the high Tc 
superconductors-as spin liquids has been discussed at 
length elsewhergj. The importance of van Hove singu- 
larities has been dramatically confirmed by recent an- 
gle resolved photoemission measurements of the doping 
dependence of the Fermi surfacecl. The qualitative dop- 
ing dependence expected theoretically for non-interacting 
electrons is sketched in Fig. 1. For overdoped samples 
(dashed curve) the Fermi surface is closed and electron- 
like. In a non-interacting model, the Fermi surface would 
grow as electrons are added, until it reached the van-Hove 
points. A Fermi surface passing through the van Hove 
points is shown as ihe solid line in Fig. 1. Instead, the 
experimental results is that as more electrons are added, 
the material develops a gap near the van Hove points. 
States in the vicinity of the van-Hove points are pushed 
away from the Fermi surface. For "underdoped" materi- 
als no states with energies near the chemical potential are 



observed near the zone edges. States are observed near 
the chemical potential only in the disconnected regions 
along the zone-diagonal shown as solid arcs in Fig. 1. 
The existence and consequences of this "non-Luttinger" 
Fermi surface require theoretical explanation. 

The subject of non-Luttinger Fermi surfaces has at- 
tracted substantial theoretical attention. The general 
approach has been to start with fermions with a large 
(Luttinger) Fermi surface and then to invoke a physi- 
cal mechanism to open a "pseudogap" which eliminates 
part or all of the Fermi surface. Three classes of mech- 
anisms have been extensively considered: (i) quasi-long- 
ranged antiferromagnetic spin fluctuations, (ii) the "d- 
wave RVB" state, and (iii) the "staggered flux phase". 
None has proven completely satisfactory; we discuss each 
in turn. 

The logic behind the antiferromagnetic spin fluctua- 
tions approach is that static antiferromagnetic order at 
wavevector Q leads to Bragg scattering at Q which may 
open a gap over all or part of the Fermi surface. Schri- 
effer and co-workers have argued that sufficiently slowly 
varying antiferromagnetic fluctuations with a sufficiently 
long correlation length may also lead, if not to a gap, at 
least to a rather strcmg suppression of the Fermi sur- 
face density of statesd. A difficulty with this picture 
is that magnetic instabilities of a wide variety of mod- 
els have been investigated; pseudogaps have only been 
found in jiiarameter regimes leading to long-ranged order 
at r = OHCI. The essential reason is that quasistatic (i.e. 
frequency much less than temperatuteVspin fluctuations 
are required for pseudogap format ionEiEI. In metallic and 
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superconducting high Tc materials, the spin fluctuations 
observed by NMILhave characteristic frequency scale of 
order T or greateru (corresponding to no long range order 
at T = 0) and are not sufficient to open a pseudogap in 
the models which have been considered. 
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Fig 1. Large square: Brillouin zone for fermions. Solid line: 
Fermi surface passing through the van Hove point. Dashed 
line: Fermi surface of overdoped material. Heavy arcs: region 
where Fermi surface was observed in photoemission experi- 
ment in underdoped BSCCO. Dashed square: inset showing 
phase space for gauge field fluctuations. Gauge field fiuctu- 
ations with momentum k in one of the four shaded regions 
are described by the conventional overdamped propagator; 
fluctuations with k outside of these regions are much larger 
because there is no part of the Fermi surface tangent to k. 

An alternative model for the origin of the pseudogap 
is the "d-wave RVB" model. This is a mathematical ex- 
pression of Anderson's original insightQ that in models in- 
volving a strong on-site repulsion and a density near one 
electron per site, singlet pairing and antiferromagnetism 
do not differ much in their short-ranged correlations and 
energies. To implement Anderson's idea one assumes 
that "spin-charge separation" occurs, and that the spin 
degrees of freedom are described by charge fermionic 
"spinous" . These fill a Fermi sea with a large (Luttinger) 
Fermi surface. This Fermi sea of spinons may be deri|Sfed. 
in a mean field theory, "the uniform RVB phase" p'ES 
of the t-J model. The t-J model is believed by manyE 
but not alH authors to contain the essential physics of 
high- Tc superconductors, and the uniform RVB phase is 
believed to be the most physically appropriate starting 
point, at least for materials near the optimal doping. To 
understand pseudogap formation one then considers in- 
stabilities of the uniform RVB state, which are due to 
residual interactions neglected in the mean field theory. 
Before now calculations have been based on "generic" 
Fermi surfaces without van Hove points and have consid- 



ered instabilities to antiferromagnetism, to d-wave pair- 
ing, and to a staggered flux phase. 

The antiferromagnetic instability has been considered 
and found not to lead to a pseudogap fer|-essentially the 
same reason as in the Fermi liquid casal3llj. 

The staggered flux phase involves the appearance of 
circulating spin currents. This entails spontaneous break- 
ing of time reversal symmetry which has been observed 
not to occur in cupratestfl. For this phenomenological 
reason the staggered flux phase has been discarded. 

The d-wave RVB state naturally leads to the forma- 
tion of a pseudogap, and the resulting phenomenology 
provides an attractive scenario for the cuprateal3. The d- 
wave RVB state may be viewed as arising from a pairing 
instability of the uniform RVB state; the resulting the- 
ory is very similar to that arising from conventional su- 
perconducting pairing. One important difference is that 
because the pairing involves chargeless "spinons" it does 
not lead to superfluidity or indeed any other observable 
which could serve as a order parameter. For this reason 
fluctuation corrections convert the pairing transition to 
a smooth crossover. 

A difficulty witii this scenario has been pointed out by 
Ubbens and Leeli3. Their results, we believe, are most 
simply interpreted as saying that the spinon-gauge-field 
interaction produces a very short inelastic lifetime for 
the spinons. This inelastic scattering is so strongly pair- 
breaking that it completely suppresses the d-wave pairing 
instability. Of course, a_£rst order transition to a paired 
state would be possiblell3, but is not observed. 

Very recently, a modeLwith an SU{2) x SU (2) symme- 
try has been consideredllZI. This symmetry is in principle 
broken down to SU{2) only at any non-zero doping, but 
the symmetry breaking is argued not to be important. 
In this model a staggered fiux instability occurs which 
does not break time reversal symmetry and can be trans- 
formed using an operation in the SU{2) x SU{2) group 
to a d-wave pairing state. The effects of gauge fluctua- 
tions on this state have not yet been determined. On the 
mean field level it leads to a phenomenology very similar 
to the one we shall derive in this paper. 

All of the previously discussed theoretical calculations 
were based on "generic" models which did not contain 
van Hove sigularities. In this communication we show 
that when the effects of proximity to van Hove singulari- 
ties are properly included, the theoretical predictions are 
very substantially altered. Most notably, near Van Hove 
points the inelastic scattering is suppressed, permitting 
the formation of a d-wave paired state with a gap which 
is non-zero only near van Hove points. This d-wave state 
should be reconsidered as a possible explanation for the 
pseudogap. 

Near a van Hove singularity the fermion density of 
states diverges, so that even arbitrarily weak interactions 
can produce large effects. This line of thinking has led to 
a large literature devoted to analyzing the implicpiions of 
van Hove singularities in Fermi liquid modelsliMcll. We 
have already mentioned the strong phenomenological evi- 
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dence for regarding high Tc materials as spin liquids^. An 
additional difficulty with Fermi liquid based treatments 
is that the singularities due to the van Hove points are 
cut off by hopping in third dimension. The bare value 
of this hopping may be determined from band struc- 
ture calculations and is not smallEj. Evidently in the 
actual materials the interplanc coupling is renormalized 
down to a very -small value in BSCCO and underdoped 
Y Ba2Cu'iO'j-:^^ this renormalization cannot be under- 
stood in a Fermi liquid picture but follows naturally in a 
spin liquid. 

The remainder of this paper is organized as follows. In 
Section II we formulate and solve a model of a spin liq- 
uid with a van Hove singularity coupled to a gauge field. 
In Section HI we employ a leading logarithm approxi- 
mation to determine the effects of residual non-singular 
interactions. Section IV is a conclusion which summa- 
rizes the approximations employed, results obtained and 
consequencies for photoemission and other physical prop- 
erties. An Appendix gives derivations of some results 
used in the body of the paper. 



II. SPINON-GAUGE MODEL 

The spin liquid state we shall start from is the "uni- 
form RVB" state. The Hamiltonian is 



t = QJ + handS 



t' 



(3) 
(4) 



Here is a number of the order of unity, in the large 
N hmit e = 2/7r2 It is believed that i' < in high 
Tc superconductorsES. Note that there is no coherent 
hopping of spinous in the third dimension; in the model 
such hopping cannot occur unless the band-structure 
ti_ ^ tband- The derivation of the spinon Ep from the 
t — J model is outlined in Appendix A. 

Near a van Hove point (e.g. px — tt, py — 0) we have 



ep = -uoipl - 
with uq = {t + 2t')a'^ and 
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If t' 7^ then 
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a 7^ 1 and the energy contour which passes through the 
van Hove points is not nested. An example of such a con- 
tour is shown as the solid line in Fig. I. At the van Hove 
points the velocity vanishes, implying a diverging den- 
sity of states and a vanishing of the fermion-gauge-field 
coupling. 

We now consider the fermion-gauge field interaction. 
Previous workEj has shown that it is correct to analyze 
this interaction in two steps using a loop expansion con- 
trolled by the parameter TV, the fermion spin degener- 
acy. In the physical problem N ^ 2. It has been shown 
that results obtained at leading order in a 1/N expan- 
sion are not significantly changed at higher orders. First, 
one constructs the renormalized gauge-field propagator 



pka^kCpj^k/2,A^p+k/2+^p-k/2]cp-k/2,a- (1) Q by drcssiug the term ff^^/gl by the fermion trans- 
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This Hamiltonian is the usual one (see, e.g. Era), however 
we have written the fermion-gauge-field coupling in its 
general form. Here v{p) = dep/dp and /^^ = tpvd^i,au- 
It represents the low energy spin degrees of freedom of 
the t — J model. Low energy in this context means en- 
ergies less than J, which in high Tc materials is about 
0.15 eV . The coefficient g^'^ contains the contributions 
to the gauge-field stiffness from the higher energy spin de- 
grees of freedom which were integrated out in the deriva- 
tion of Eq. (|l|) . It may be expressed in terms of the high 
energy part of a six spin correlation function and is of 
the order of J in magnitude; it is discussed further in 
Appendix A. 

We assume the fermion spectrum is 



-2t[cos(pxa) + cos(pya)] 

-At' cos,{pxa) cos{pya) — ji — At' 



verse current-current polarizibility H = J vvGG. Sec- 
ond, one uses this to compute the fermion self-energy. 
For a closed Fermi surface far from van Hove points, 
■jT{uj,k) ~ ^^+xfc^- For a given direction of k, the dissi- 
pative term |w|/|fc| comes from fermions near the points 
on the Fermi surface which are tangent to k; po is the 
curvature of the Fermi surface at these points and x is 
the diamagnetic susceptibility of the fermions. 

The presence of van Hove singularities leads to two ef- 
fects. First, as can be seen from Fig. 1, there is a range 
of directions of k which are not tangent to any point on 
the Fermi surface. For these directions, the dissipative 
term in Il(LU,k), Ildiss, is much smaller. To display the 
behavior of the dissipative term it is convenient to define 
coordinates k± = k^ aky . 

Ifl^<|A:+A„|<i±^,then 



(2) 



aN\uj\ 



Here a is the lattice constant, fi + At' is the chemical po- 
tential with /i = corresponding to the van Hove point, 
i is a first neighbor hopping due mostly to the superex- 
change J with an additional contribution coming from a 
bandstructure tband renormalized by the hole density and 
the parameter t' is derived from further neighbor hopping 
in the underlying bandstructure also renormalized by the 
hole density; roughly, we expect 



fc2 



fc2 



Att \k^ ~ a^k^ k^ — a^k"^ 



while if is outside this range then 



Npo\uj\ 
27r|fc| 



(6) 



(7) 



Second, states in the vicinity of the van Hove point pro- 
duce a negative, divergent contribution to x, so 
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Fermi surface; from Eq. (^) we conclude that 1/po oc 5, 
so CJq oc S. 



Here the logarithm cutoff A = max(T, fj,,uok^) and 
Xreg IS the Contribution from fermions far from the van 
Hove points. The total x calculated for non-interacting 
fermions with N = 2 moving in the band structure de- 
fined by Eq. (j^) is shown in Fig. 2; as seen in this Figure 
X near half filling is negative and of small magnitude with 
a weak divergence near van Hove points as expected from 
Eq. (^). The full gauge propagator is thus 



D{uj,k) 



1 



+ (.9o +X)fc' 



(9) 



The divergence of x implies that in the vicinity of the 
van Hove singularity the uniform RVB state is unstable 
to a state of non-zero flux, at a temperature 



flux ~ exp ^ 



Nuo 



(10) 
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unless a different instability occurs first. States of non- 
zero flux break time reversal invariance; as there is no 
evidence that this occurs in high Tc materials we shall 
assume that Xreg + ffo"^ is sufficiently large that Tfiux is 
negligibly small. This assumption is consistent with the 
previously mentioned theoretical estimate ^ J which 
is greater than the typical values of x ~ 0-1^ ^ O.IJ 
shown in Fig. 2. Also, a previous analysis of the T- 
dependence of the resistivity at T > 100 K predicted by 
Eq. (|l|) found g^^ + Xreg « 500a2 K - These es- 

timates combined with the large numerical factor in Eq. 
( p^ imply that the instability will occur only very close 
to the van Hove point, and only at very low temperature. 

At scales of interest we may then neglect the logarithm 
and write 



D{uj,k) 



1 



(11) 



where g~^ —90 + X ~ 500a^ K. We now use this 
D{uj, k) to calculate the fermion self energy. For fermions 
far from the van Hove points the, calculation is identical 
to that given in previous worknJ. The kinematics of a 
scattering event imply that the fermion is scattered par- 
allel to the Fermi surface. For these momenta one must 
use the |tf form of Xldiss, leading to 



1/3^2/3 



Here 



(2V3)3 



7r2po 



(12) 



(13) 



at leading order in N . The Fermi velocity vpip) vanishes 
linearly as one approaches the van Hove point, implying 
S far does also. Note, also that 1 /po vanishes for a nested 



Fig 2. Diamagnetic susceptibility of non-interacting Fermi 
gas with spectrum given by Eq. (2) as a function of filling, p, 
for t = 1 and several values of the next neighbour hopping, 
t'. 

For fermions near the van Hove point, scattering pro- 
cesses involving the other form of Hdiss become allowed. 
These lead to a contribution 



In — In ■ 



[iAe-Bu{pl 



1 

27r^ V e Ti'Xy 

(14) 

where the numerical coefficients A = In — 1, B = 
Snear(e,p) bccomes important for e < e* 
In — K, 1, so In ^ is within 
thus, for consistency, we 



31ni±ii - 1 

1 — a 



where -J- 



In^ln^ = 1: 
a factor of two of In — ^ 



must assume I]„ear(e,p) is negligible. The instability in 
S„ear(e,p) reiuforccs that in x- To study this instability 
more carefully we have written and solved coupled RG 
equations for u and g; we find that the more careful con- 
siderations do not change the estimate of the scale Tfi^x 
where these effects become important. Thus we assume 



S(e,p) = S/ar(e,p); 



(15) 



in other words the self energy due to fermion-gauge field 
scattering is important far from the van Hove points and 
unimportant near them, provided that the instability to- 
wards the flux phase may be neglected. 



III. EFFECT OF NON-SINGULAR 
INTERACTIONS 

We now turn to the residual, non-singular interactions. 
It has already been shown that the short range interac- 
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tion between fermions far from V|aiij-IJove points is renor- 
malized to zero by the gauge fielcCjO. We thus need only 
to consider processes involving fermions near van Hove 
points. For these processes, vertex renormalizations due 
to the gauge field lead to the same logarithms we have 
already agreed to neglect. The theory we must consider 
is therefore of fermions near the van Hove points, with 
self energy I]/ar(e,p), coupled by short range interactions 
which may be parametrized by four constants: particle- 
particle and particle-hole interactions with momentum 
transfer near zero or G = (tt, tt). 

The theory has some formal similarity to theories of 
fermions in one dimension. In one dimensional models 
one considers only fermions near the Fermi surface of the 
left and right hand branches of the dispersion relation; 
here one considers only fermions near van Hove points. 
In one dimensional models a weak coupling leading log- 
arithm approximation exists because some particIe-hoIe 
and particle-particle propagators diverge logarithmically 
due to kinematics in one dimension. In the present model 
the divergent density of states at the van Hove singularity 
similarly leads to logarithmic divergences of susceptibili- 
ties. 

The theory is more complicated than theories of 
fermions in one dimension because there are two small 
scales: Wq = Tr^ep which will be shown below to set the 
scale at which deviations from the perfect nesting become 
important, and luq which sets the scale at which gauge 
field effects become important. Both cua and loq are pro- 
portional to 6 as previously discussed; because they have 
the same doping dependence and ujq is small even at large 
doping due to the large value of the gauge stiffness 
in Eq. (^3|), we believe that ujq < uja is the only relevant 
case. 

The charges in the antiferromagnetic and rf-wave pair- 
ing channels diverge. At scales larger than LUa and loq 
both channels have a ln^(l/e) divergence. At scales 
less than Ua the antiferromagnetic divergence becomes 
In^l/e). At scales less than ojo the coefficient of the 
In (1/e) in the superconductivity channel (and, if uq > 
LOa, in the antiferromagnetic channel) is reduced, because 
the gauge field produces a strong inelastic scattering in 
some regions of momentum space. 

In the remainder of this Section we present the re- 
sults of the leading logarithm calculations. We assume 
throughout that at most one coupling becomes large. Es- 
pecially in the regime of ln^(l/e) renormalizations the 
problem of coupled, charges is very involved and has been 
treated elsewhera22l. 



ID, the small momentum transfer particle hole suscep- 
tibility has only a weak singularity as w = vq, so these 
processes can be neglected. The G momentum transfer 
particle-hole processes lead to logarithmic divergences. 
For a nested Fermi surface (a^ — 1) the two momen- 
tum integrals are each logarithmically divergent leading 
to a ln^(l/e) renormalization of the charge. For a non- 
nested Fermi surface (a^ 7^ 1) one logarithm is cut off by 
(1 — a^) if it has not already been cut off by temperature. 
The crossover occurs at the scale w « . The existence 
of a logarithmic divergence in the particle hole bubble 
means that the charge, gAF, associated with G momen- 
tum transfer particle-hole processes grows. In the leading 
logarithm approximation we find that antiferromagnetic 
charge is renormalized by 



S9AF _ r m: 



OAF 



9AF 1„2 ep_ 
gAF l+a ep^Q 



2^2 



(16) 



The second of these formulas was derived on the assump- 
tion that u>a > Wo; if not an additional crossover occurs. 
The formulas relevant to this case will be presented and 
discussed in the next subsection treating pairing. The 
renormalization of qaf calculated from Eq. (^|) is shown 
as the dashed line in Fig. 3. 
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Fig 3. Renormalization of the interaction constants g = 
^ plotted against t = \n{eF/T). Sohd hne: ^ calcu- 

9p 



9AF ' Bp 

lated from Eq. (I 



with initial condition 



3p 

0.1. Dashed 



line: 



9AF 

Mf- = 0.2, 



calculated from Eq. ([L6|) with initial condition 



A. Antiferromagnetic instability 

We proceed to construct leading logarithm RG equa- 
tions. There is no logarithm in the small momentum 
transfer particle-hole susceptibility, so the correspond- 
ing charge is not renormalized. Moreover, in contrast to 



We define Taf as the scale at which -2ae. becomes of 

gAF 

the order of unity. At scales larger than Taf , the effects 
of gAF are negligible within the weak coupling approx- 
imation. At scales of order Taf, the antiferromagnetic 
susceptibility becomes large, the renormalization of the 
fermion propagator becomes large, and perturbation the- 
ory breaks down. It is difficult to make definite predic- 
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tions for temperatures lower than Tap because one is 
then deahng with a strongly interacting model. We can 
imagine two scenarios - a crossover to the d-wave RVB 
regime to be discussed below or a crossover to an anti- 
ferromagnetic regime involving critical fluctuations. We 
discuss the antiferromagnetic regime here, focussing on 
whether the important fluctuations are quantum or clas- 
sical and whether they can open a gap in the fermion 
spectrum. 

It is helpful to compare the present calculation to the 
well known BCS theory of superconductivity which also 
has a logarithmic divergence of the coupling gsc- In su- 
perconductivity the important fluctuations are classical 
(i.e. involve modes with energies less than ksT) and are 
generally weak but grow as T approaches Tc- These fluc- 
tuations are described by the classical Landau theory 



(17) 



where ly ^ l/ep is the density of states, the coherence 
length ^0 = vp/Tc and f^c- The renormalization 

of (3sc may be estimated to beEa 



5^ 



(18) 



where r 



This estimate shows that at the scale 



T « 1 where the coupling g^c has scaled to the order 
of unity, thermal fluctuations are negligible in dimension 
d > 1. By continuity quantum fluctuations must be also 
negligible at this scale. This may be seen directly: at 
scales larger than Tc, the cj)'^ vertex of the quantum Lan- 
dau theory /3sc{e, k) ~ e^+(\>fc)^ ^^'^ ^^"^ Cooper propaga- 
tor has negligible momentum dependence so the leading 
correction SPsc/Psc ^ (Tc/^f)'^^^- Thus we see from 
the quantum calculation that d = 1 is the marginal di- 
mensionality for the quantum fluctuations; this fact is 
revealed in the classical calculation because in d = 1 the 
Ginzburg parameter (Tc/ep)'^'^ ceases to be small. 

In the antiferromagnetic problem of interest here, the 
marginal dimension is d = 2 because the fermion spec- 
trum is e oc rather than e oc p. The free energy for the 
classical fluctuations is 



In 1+^ 



Fit 



'Tai 



'AP 



(19) 



with (3af = Fs/iuTlp) and r = {T- Tap) /Tap {F, are 
numerical coefficients of the order of unity). From the 
free energy ([l9|) one finds 



AF 



f3AF 



G,. = In- 



1 



a 



1 



(20) 



We see that d = 2 is the critical dimension because the 
only parameter (apart from t) controlling the fluctuation 



correction is In which does not depend on Tc/ep- 
The calculation of the leading quantum fluctuation cor- 
rections shows that these also are small only by a power of 
In-i^^. We cannot proceed further in the general case, 
but in the limit a ^ 1 we can use the In to con- 

1 — Q 

trol the calculation assuming that the antiferromagnetic 
charge is the only relevant one. However, a theory of the 
general a ^ 1 case found more than one relevant charge 
and the resulting theory of the competing instabilities is 
very involvedEll. 

In our a — )■ 1 limit, the only important fluctuations 
are classical. There is no long range order at any T > 0; 
rather, for T < Tap, A crosses over to A ~ exp (^) im- 
plying an exponentially growing correlation length ^ap ^ 
^-1/2^ These fluctuations are thus quasistatic and long 
ranged and in particular have energy much less than the 
typical fermion energy fegT and momentum much less 
than the typical fermion momentum, p ^ ^JT ju, so as 
far as the fermions are concerned these fluctuations may 
be treated as static periodic scatterers, and lead to a 
gap Ai? « Tap^/ —T . The magnetically induced fermion 
pseudogap found in this model is not relevant to the un- 
derdoped high Tc materials because the fluctuations pro- 
ducing the pseudogap would also lead to a very rapid 
T-dependence of the Cu NMR relaxation rates propor- 
tional to powers of £^ap ^ 'sx^Y'iGiTAP /T); such a rapid 
temperature dependence is not consistent wittuCu relax- 
ation rate measurements in high Tc materialsa. 



B. Superconducting instability 

We now consider particle-particle processes. As usual 
the leading divergence happens in the q — Q momentum 
channel. The Cooper propagator, C{T), is 



CiT) 



(Ppde 



9p 



(27r)3 e2 + (g + ^^(p)l/3£2/3)2 



(21) 



Near the corners ujo{p±) 
so one finds 



[p+p-f and ep = -up+p-, 



C{T)', 



Sp 1 2 tf. 
8^ ^ T 

247r^M I 



In^ 



-61n ^ 



T >ujo 
T <ujo 



(22) 



Here gp is the pairing coupling constant. 



The In ^ divergence of C (T) is due to the divergence 
of the density of states and the vanishing of the inelastic 
lifetime near the corners. The change in C(T) as T is 
reduced below coo reflects the pairbreaking effect of the 
gauge fluctuations in the regions away from van Hove 
points. One sees immediately from (|2l| ) that C(T) is 
non-divergent in models where uJo{p) ~ const, as shown 
by Ubbens and LeeEEl. Ubbens and Lee derived their re- 
sults by an energy argument which demonstrated that a 
second order pairing transition is not possible in models 
with ll>o{p) const. We believe the argument using the 
gap equation is equivalent but more transparent, however 
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we note that if we apply their energy arguments to the 
present model we discover that a second order transition 
is possible. 

The In^ ^ divergence of C{T) means that pairing is 
the instability which dominates the logarithmic antifer- 
romagnetic and flux instabilities in weak coupling. How- 
ever, because the interactions are repulsive the pairing 
instability occurs in the d-wave channel. Formally, one 
must consider the pairing amplitudes An and A22 near 
the two inequivalent corners 1 and 2 and the appropri- 
ate interaction amplitudes gn, gi2. The gap equation 
is dominated by the corners, so becomes a 2 x 2 matrix 
equation 



An 
A22 



-C{T) 



911 512 
512 511 



An 
A22 



(23) 



The growing eigenvalue is for An — — A22 (corre- 
sponding to d-wave symmetry) and requires that gp = 
.912 — gii > 0. The pairing scale Tp at which the interac- 
tion becomes of order unity is given by 



e F exp - 
4 

exp- 



5p 



61n- 



2 



Tp > loq 
Tp < luq 



(24) 



In a model with point-like interactions, 1712 = 511; it 
is then necessary to go to higher order. Contributions to 
the pairing instability of order 1ii'^(c-f/^) exist and, 
for sufficiently small coupling dominate over the AF in- 
stability. The renormalization of gp is shown in Fig. 3. 

By repeating the arguments of the previous Section we 
see that the Ginsburg parameter is Gi ~ l/ln{eF/Tp), so 
mean field theory provides a reasonable description at 
temperatures less than Tp in weak coupling. There is, 
however, one important caveat: mean field theory pre- 
dicts a phase transition to a paired state Tp] in reality 
this is not a true transition. The interaction between 
vortices of the pairing field is not logarithmic because of 
screening by the gauge field, so as long as spin-charge 
separation occurs vortices exist and prevent long range 
order. However, the small value of the Ginzburg param- 
eter implies that the number of vortices present is very 
small, so the magnitude of the pairing field and the im- 
plications of the pairing for the physical properties are 
well described by the mean field theory. In particular, in 
regions of momentum space where the pairing amplitude 
is appreciable there will be a strong suppression of the 
electron spectral function. 

We now discuss the momentum space structure of the 
gap function focussing on the physically relevant Tp < coq 
case. The gap equation is 



gpair{i,p)At',p'{dp') 



ie' + 4^'{p')ey^r + {ep,+A,,yy 



(25) 



The integral in ( p5| ) is dominated by the van Hove point 
which allows us to ignore the p\ e' dependence of gpair in 



(p5|). For p away from the van Hove points .gp(e,p,p'Vis 
suppressed at low scales by the gauge field fluctuationsEJ, 
roughly gpair{e,p) = Opizj^T for e < wo(p) with 
K > 2/3. Eq. (1^) implies 



A 



wo(p) < A 
ujq{p) > A 



(26) 



Here A ~ Tp is the pairing amplitude at the van Hove 
points. Clearly, as one moves far enough from the van 
Hove points so that tooip) > A, the pairing amplitude 
drops rapidly and becomes less than the scattering rate 
ujQ^^{p')e'^^'^ due to the gauge field; these portions of the 
Fermi surface may be regarded as gapless. 



IV. CONCLUSION 

We have studied theoretically a model of a spin liquid 
with a Fermi surface which passes near a van Hove singu- 
larity. We considered two sorts of interactions: the singu- 
lar gauge interaction arising from the spin-charge sepa- 
ration which established the spin liquid in the first place, 
and residual short range interactions between spinous. 
The gauge field interactions lead to two effects. One is 
an instability ai T — Tfi^x to a "flux phase" in which 
time reversal symmetry is spontaneously broken. An ex- 
pression for Tfiux is given in Eq. (|lO|). The numerical 
factors are such that Tfi^x is negligible in the weak cou- 
pling limit. As there is no experimental evidence for time 
reversal symmetry breaking in high Tc materials, we as- 
sume that parameters in the physical model are such that 
Tfiux is negligibly low. This assumption allows us to ne- 
glect also the logarithmic renormalization of the fermion 
dispersion shown in Eq. (^5|). The second effect is a sup- 
pression, for spinons near the van Hove point, of inelastic 
scattering due to the gauge field. The physics is simple: 
the gauge field couples to the fermions via the velocity; 
this vanishes at the van Hove point and the vanishing 
coupling overcomes the diverging density of states. 

We then turned to the non-singular interactions. Pro- 
cesses involving spinons far from the van Hove points are 
known to be renormalized to zero by the gauge-field, but 
near the van Hove points the vanishing of the coupling 
implies that the renormalization is ineffective. We there- 
fore argued that we could specialize to a model involving 
fermions near van Hove points coupled by non-singular 
interactions. We treated this theory via weak coupling 
leading logarithm methods similar to those used to study 
one-dimensional models. 

We found diverging interactions in the d-wave pairing 
and antiferromagnetic channels. The coupling constant 
flows are shown in Fig. 3. For sufficiently weak cou- 
plings, d-wave pairing dominates and a controlled ex- 
pansion based on the parameter ln(ei?/T) is possible. 
Below the pairing temperature Tp given in Eq. ( |2^ ) 
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the d-wave pairing leads to a gap in the fermion spec- 
trum near the van Hove points but leaves a finite re- 
gion of gapless Fermi surface near the zone diagonal. 
The resulting fermion spectrum is similar|-|to that ob- 
served in recent photoemission experimentsu. This pair- 
ing which eliminates some but not all of the Fermi sur- 
face will have implications for other physical properties. 
For example, the uniform spin susceptibility will decrease 
as T is decreased through the pairing scale, but will 
tend to a non-zero limit as T — > because some of 
the Fermi surface remains ungapped. Precisely this be- 
havior occurs in La2-xSrxCuQi at an x-dependent pair- 
ing scale T*{x) varying from T*{x = .15) « 300 K to 
T*{x = .04) « 700 Other properties also exhibit 
crossover at T*(x)tZI and similar behavior may pecur in 
Y BaCu^O^+x^ although the data are less cleared. The 
identification of the RVB pairing scale Tp with thCjeHipir- 
ical T*{x) was first suggested by Tanamoto et a£3; the 
new features of the present paper are a theoretical justifi- 
cation for the existence of the gap despite strong inelastic 
scattering and the result that the gap opens only over a 
small portion of the Fermi surface. 

The present theory does not explain the rapid drop 
of Xs{T) and NMR relaxation rates observed in bilayer 
and trilayer materials below a spin gap temperature 
TsG ^ 200-?^. Explaining these observations requires 
opening a gap over the whole Fermi surface; Oiie pos- 
sibile mechanism has been discussed elsewherecj. Fi- 
nally, we note that in this theory the d-wave RVB state, 
and therefore the fermion pseudogap, appears at a dop- 
ing 6g > 6vH- Experimentally the gap appears first for 
Sg ~ 6vH and it is important to determine if Sg > Svh 
or not. 

The theoretical situation is less clear if the antiferro- 
magnetic channel is dominant. There is no generic small 
parameter to control fluctuations in this regime. How- 
ever, if the antiferromagnetic channel is dominant and 
the Fermi surface is nearly nested, a controlled calcula- 
tion turned out to be possible. In this case we showed 
that the physics is controlled by classical spin fluctua- 
tions with correlations which grow exponentially as T 
decreases below a mean field scale Taf- These qua- 
sistatic fluctuations produce a gap, A^?, in the fermion 
spectrum near the van Hove points which varies as A^? = 

Taf \J '^'ta~'^ ' foi' ^ < ■ This mechanism for produc- 
ing a pseudogap predicts a very rapid T-dependence of 
NMR rates which is not- observed in experiment on su- 
perconducting materialsU. 

We now discuss the doping dependence of our results. 
In a rigid band model one would expect Tp to rise and 
then fall as the chemical potential is tuned up to the 
van Hove point and then beyond. However, in the spin 
model the particle-hole symmetry breaking term t' itself 
depends on doping as seen in Eq. (^. As the doping 
is decreased beyond the van Hove point the Fermi sur- 
face gets flatter. This tends to pin the Fermi level to 
the corners. Further, as the Fermi surface flattens, po in- 



creases, so that Wo in Eq. ( |T^ decreases and Tp increases. 
However, Ua also decreases so the importance of antifer- 
romagnetism grows. As the doping is reduced towards 
zero, all instability scales become greater than both ujq 
and uja ■ In this regime the pairing and antiferromagnetic 
charges scale in the same way, as shown in Fig. 3. Be- 
cause we expect the bare value of the antiferromagnetic 
charge to be larger than the bare value of the pairing 
charge we expect antiferromagnetism to be dominant, as 
was also found by Dzyaloshinshkii and Yakovenkocil. 

Within our leading logarithm approximation we have 
shown that the observed pseudogap can not be due to 
long-ranged antiferromagnetic fluctuations, and must be 
due to RVB pairing. We now offer qualitative arguments 
that this conclusion survives even the leading logarithm 
approximation is not reasonable. First, it seems that 
antiferromagnetic fluctuations can produce a pseudogap 
only if they are long ranged and quasistatic; such fluctu- 
ations are ruled out by NMR so we believe that a pairing 
origin of a pseudogap is more likely. Second, increases in 
the antiferromagnetic charge seem to feed back into the 
pairing equation in a way that increases Tc, so it seems 
natural to expect that the theory with antiferromagnetic 
charge ~ 1 is unstable to d-wave pairing if the Fermi 
surface is not too flat. However, when the Fermi surface 
becomes flat, antiferromagnetism becomes favored. On 
the basis of these considerations we propose the phase 
diagram shown in Fig. 4. 
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Fig 4. Phase diagram of the spin liquid near van Hove 
point. AF denotes antiferromagnetism, SL denotes spin liquid 
with ungapped Fermi surface and PG denotes the pseudogap 
regime of the spin liquid phase 

Finally, we discuss the spin fluctuations expected in 
a different regime of this phase diagram. In the large 
(5, spin liquid, regime and in the pseudogap regime near 
the SL-PG boundary we expect spin fluctuations-|dom- 
inated by the "2^^" effects discussed elsewherai3. As 
S is decreased we expect enhancement of the spin fluc- 
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tuations near (tt, tt) point due to the flattening of the 
Fermi surface and the increase of the renormahzed inter- 
action. Moreover, inside the PG regime there will be a 
gap on those parts of the Fermi surface which could damp 
a (tt, tt) fluctuation. Therefore these fluctuations would 
be undamped as is apparentlv-^equired by analyses of Cu 
NMR Ti and T2 experimental. 

In summary, we have studied the effect of van Hove 
points on the "uniform RVB" state with weak residual 
interactions. We have identifled three instabilities: to 
a flux phase with spontaneous breaking of time reversal 
symmetry, magnetic phase, and a d-wave RVB phase. We 
argued that the d-wave state is likely to happen at large 
doping while antiferromagnetism dominates at small dop- 
ing. Our proposed phase diagram is shown in Fig. 3. We 
showed that the inelastic scattering due to the gauge fleld 
is negligible for fermions near the corners and that this in 
combination with the divergent density of states permits 
a continuous crossover to a d-wave RVB state, which is 
not allowed in models without corners. 



We now estimate the bare gauge stiffness go . This rep- 
resents the effect of short wavelength fluctuations which 
are integrated out in the definition of Eq. (^; in the 
N 00 limit gQ^ vanisheso, but it is non-zero for the 
physical N — 2. If the only term in H, Eq. (|^), were 
f^^, then one would find 



1 

4sJ 



dt{h,{t)h,{0)) 



9l 



(A5) 



Here as usual h = dxCLv 



To find (7o we evaluate the 



{hi(t)hi{0)) correlator directly from S = 1/2 Heisenberg 
model, including in our calculation only short wavelength 
degrees of freedom. Formally, the gauge field is related to 
the phase of the operator Aj_; fWhich decouples the four 
fermion operator in Eq. ( |A2| )d; the field h is equal to 
the gauge flux through the elementary placquette with 
vertices i,j,k,l; the latter is related to the A operators 

by 



(A6) 



APPENDIX A: PARAMETERS OF THE SPINON 
HAMILTONIAN 

The spinon Hamiltonian is obtained expanding about 
a mean fleld approximation to the t— J model which may 
be written 



Here i? is a real operator. By undoing the decoupling 
one may express the product of four A around the plac- 
quette in terms of fermion operators; in the low doping 
limit these in turn may be expressed in terms of spin 
operators (with corrections 0((5)), yielding 

+ ^l:i^,k)S^{S,y.Su) (A7) 
+2[{S,Sj){SkSi)+{Sji) {SkS, )~iS,Sk) [Sj Si )] 



Here is the spinon creation operator discussed in 
the text, Si = ^J2ai3 4a^ai3Cia, thc bosc operator 6| 
creates a spinless charge 1 "holon" and the constraint 
^l^i + J2a ^a'^ia = 1 is assumed. To obtain the fermion 
dispersion one writes the second term in Htj as a product 
of four fermion operatoraa 

SiSj — ^-^CjaCiacljjCjp + - (A2) 

and approximates this by 

SiSj « -^(c]„c„ + h.c.){cl^Cjp). (A3) 

If {c\i^Cjp) ^ then the bosons acquire a dispersion 

Hhose — ^{c\pCji3)t\°:^'^h\hj] because the density of 
bosons is low the bosons are mostly in small k states 
so it follows that for nearby i and j {b\bj) — 6. Inserting 
this in Eq. (Al) we get the fermionic part of Eq. (|l|) 
with 



{band) 



(A4) 



The expectation value {clpCjp) turns out to be -^+0{6). 



where (ij) denotes all distinct pairs of spins, (ijk) de- 
notes all distinct triads with the clockwise order of k 
around the placquette. If there are strong short-range 
antiferromagnetic correlations then we may replace the 
real terms by c-numbers, obtaining 



Re 



ih 



1 



E S^{S, X Sk) 



This allows us to identify 



h^2j2 S^iSJ X Sk) 
(ijk) 



(A8) 



(A9) 



and thus to relate g§ to a six spin correlator. Evaluating 
this in the spin wave approximation leads to 



9l 



1 



k,p 



tOk 



(AlO) 



where i? is a number of the order of unity and and 
LUp are spin wave energies. Because the typical value for 
LUk ~ 2J and the sums are dominated by short wave- 
lengths, we see that g^^ = B'J with B' a number of the 
order of unity. Of course we can not expect this simple 
estimate to yield a reliable value for B', but it does show 
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that it is reasonable to expect a large contribution to the 
gauge stiffness from short range correlations. 

Finally we consider the effective between planes hop- 
ping t±. If ^ then also the b etwe en planes exchange 
J_L ^ and we must add to Eq. ( |AlD the terms 



-Y.^j^''^^'^ +^T.^rsr (All) 



•(1)^(2) 



One may factorize the four fermion term as in Eqs. 



(A2,A3); however the equation determining the ampli- 
tude is 



{ie-tpf 



(A12) 



A solution only becomes possible for J± + 6t± « t, so 
for physically relevant parameters there is no coherent 
between planes hopping in the model. 
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